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factor one. The chemical potential of such a state scales with the density as /i oc n 3 ^ 2 . This is a 
lower energy per particle than n oc n for the earlier suggested groundstate candidates: a condensate 
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I. INTRODUCTION 

Precise control over interactions of ultracold atoms 
with laser fields opened an opportunity of fabricating 
synthetic non-Abelian gauge fields 1 - - — and spin-orbit (SO) 
couplings^ 5 - of both Rashba£ and Dresselhaus^ type. In 
analogy with previously known vector condensates^—, 
an important tool for synthesizing such systems is the 
controlled Raman coupling between oriented sequence 
of initially degenerate hypcrfme states of 87 Rb. Inter- 
actions of spatially varying laser field with hyperfine 
atomic levels effectively produces synthetic non-Abelian 
gauge fields ; 13 ' 16 " — originating from the Berry phase, and 
SO couplin g 4 ' 5 ' 14 ' 15 ' 1 9 ~ — of momentum to the internal 
isospin degrees of freedom. Bosons with SO coupling, dis- 
cussed and studied in the context of cold-atom physics 
in Ref. [l5j . were realized very recently by Spielman's 
group at NIST— . Depending on the particular experimen- 
tal scheme, i.e. the choice of atomic states and a sequence 
of optically induced transitions between them, one may 
in principle realize various SO Hamiltonians for the pro- 
jected low-energy states. Low energy properties of such 
SO bosonic systems have been discussed in Refs. fl5l. [23l— 

n. 

While fcrmions with SO coupling were extensively 
studied over the last decades, see e.g. Refs. [3^ - l37| for 
reviews, relatively little is known about SO bosons. Yet, 
they offer a number of fundamental problems which, in 
some respects, are more challenging than their fermionic 
counterparts. Indeed, in many instances SO coupling 
leads to single-particle dispersion relations which exhibit 
multiple minima or even degenerate manifold of mini- 
mal energy states. Fermionic many-body systems form a 
unique Fermi sea state on top of such dispersion relation. 
On the other hand, bosons tend to occupy the lowest en- 
ergy states and thus face macroscopic degeneracy of their 
non-interacting groundstate. It is entirely the effect of 
collisions (i.e. boson-boson interactions) which lifts this 
degeneracy and selects a true many-body groundstate. 

In this respect the problem is somewhat similar to frac- 
tional quantum Hall effect (FQHE). There the fermionic 
kinetic energy is degenerate due to Landau quantiza- 
tion and the nature of the groundstate is determined by 



the interactions. One of the most remarkable concepts, 
which emerged in FQHE studies, is that of composite 
particles 4 ^— . For example, FQHE states with filling 
fractions v = p/(2p + 1), where p = 1,2..., were under- 
stood as integer v = p quantum Hall states of composite 
fermion particles. The latter are obtained by binding the 
original fermions with flux tubes carrying exactly two flux 
quanta. Technically such a binding is achieved by assign- 
ing Chern-Simons (CS) phase factor to the many-body 
wave function of composite particles. 

The goal of this paper is to show that a very similar 
transformation plays a major role in understanding of the 
groundstate of 2D bosons with Rashba SO coupling. We 
argue here that their groundstate wavefunction may be 
approximated by that of the Fermi gas at integer filling 
factor v — 1, dressed by the Chern-Simons phase with 
one flux quantum attached to every fermion. Such phase 
factor transforms the integer quantum Hall fermionic 
wavefunction into a bosonic one. Fermionization of the 
bosonic system allows the latter to minimize its interac- 
tion energy. This is due to the fact that fermions with 
the same spin can't be at the same spatial point and thus 
cannot interact through a short-range s-wavc interaction. 
(Consequently the amplitude for two fermions with al- 
most parallel spins to be at the same point is small as 
the angle between their spins.) As the result the interac- 
tion energy per particle in a Fermi sea is smaller than in 
a Bosc-condcnsate of the same density. For low enough 
density such reduction of the interaction energy wins over 
the associate increase of the kinetic energy. 

A very similar physics is behind the so-called Tonks- 
Girardeau limit of spinless ID Bose gas^r— . At a small 
density (which in ID is the same as strong interactions) 
the bosonic wavefunction approaches symmetrized wave- 
function of non-interacting Fermi gas. If spin degree of 
freedom is present in ID, the groundstate is known to 
be fully spin-polarized 4 ^— again leading to the Tonks- 
Girardeau construction in low density limit. We thus 
notice that 2D Rashba bosons share features of both 2D 
quantum Hall systems as well as ID Bose gases. The 
deep connection with the latter is due to the fact that 
single-particle density of states for particles with Rashba 
SO coupling behaves as e -1 / 2 at small energy. This is 
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typical for ID systems, making particles with Rashba 
SO coupling "ID-like", irrespective of their actual spa- 
tial dimensionality. 

Technically blending FQHE and Tonks- Girardeau 
ideas with Rashba SO coupling presents one with a num- 
ber of challenges. Indeed, the standard way of intro- 
ducing Chern-Simons transformatio n 45 ' 47 essentially re- 
lies on the spinless nature of particles. Here we suggest 
a way to generalize it for spinor particles with strong SO 
coupling. It achieves the goal of fermionization of bosons, 
but fails to eliminate completely the interactions between 
the composite spinless fermions. Still the fermionic in- 
teraction energy can be parametrically smaller than the 
bosonic one. The crucial observation here is that the 
residual fcrmion-fcrmion interactions appear to be pro- 
portional to the angle between momenta of two scattering 
particles. Therefore, confining the Fermi sea of compos- 
ite fermions to a small fraction of the momentum space, 
allows one to lower their interaction energy. Similar ideas 
were recently put forward by Berg, Rudner and Kivelson 
(BRK)^ 4 - in the context of Fermi gas with Rashba SO cou- 
pling. They called the resulting time-reversal symmetry 
broken state a nematic. Here we adopt their construction 
for the composite fermions. 

The paper is organized as follows: in section |H1 we in- 
troduce Rashba SO coupling on a single-particle level. 
We also discuss earlier mean-field theories for bosonic 
many-body groundstates along with BRK nematic state 
for many-body fermionic system. In section UHl we intro- 
duce fermionization of spinor bosons and evaluate the en- 
ergy of the resulting composite fcrmion state. In section 
IIVI we discuss the results: present an emerging ground- 
state phase diagram of Rashba bosons, list possible ex- 
perimental consequences and generalizations of our ap- 
proach. Some technical details of the calculations are 
relegated to the Appendix. 



II. RASHBA SPIN-ORBIT COUPLING 

A. Single particle spectrum 

The single-particle Hamiltonian with the Rashba spin- 
orbit coupling in two dimensions takes the form 

ffo = -^+n-z.[<TxV r ], (1) 
Ira 

where v is spin-orbit coupling constant having dimen- 
sionality of velocity, V r = {d x ,d y ) and er = (a x ,a y ) is 
vector of Pauli matrices acting on two component spinor 
ip(r) = (ip(r, t), ip(r, I)). In 2D the Rashba term in 
Eq. ([1]) may be transformed to another form ivcr ■ V r 
by 7r/2 rotation in the spin space. 

In addition to translational and rotational symmetries 
the Hamiltonian H$ commutes with the two discrete Zi 
symmetry operations: time-reversal T and 2D parity P. 



They may be defined in the following way 

f /^(r,t)\ ( i>{r,i)\ m 
{iir,l) ) - V-<Kr,t) ) (Z) 

and 

where bar stands for complex conjugated and we intro- 
duced the complex 2D coordinate as z = x + iy. Notice 
that the parity operation in 2D is defined as y — > — y 
and x — > x and o~ y multiplication. Indeed, reflection of 
both coordinates x and y is equivalent to 7r rotation. Of 
course, one could as well define parity as x reflection and 
o~ x multiplication. It is equivalent to the P operation J3]) 
combined with the ir rotation. 

Diagonalization of the Hamiltonian ([T]) yields the fol- 
lowing single-particle spectrum: 

£k , 7 = |-+ 7W fc. (4) 

Here index 7 = ±1 labels the two brunches of the 
spectrum depicted in Fig. Q] The corresponding eigen- 
functions are plane waves in coordinate space multiply- 
ing coordinate-independent spinor, whose form depends, 
however, on the momentum k as 

1 ( — Wp- iar s( k ) \ -, 
^ l7 (r >S ) = -^ }e- ^ 

where arg(k) is the angle between the momentum vec- 
tor k and the fc^-axis, V is the system's volume. The 
spin is directed in the x — y plane and is rotated rel- 
ative to the momentum direction by 777/2. Hereafter 
we consider SO energy scale eo = mv 2 /2 as the largest 
energy in the problem. One thus expects the ground- 
state of an interacting system to be confined to the 
part of the Hilbert space projected on the lower branch 
7 = — 1. Clearly the two spin components of all such 
states are connected by the following unitary transfor- 
mation ^k,-(r, t) = '4e~' iarg ( k )'0i (; ,-(r, \). For a generic 
wavefunction belonging to the lower branch the relation 
between its up and down components acquires the form 

tf(r,t) = J dv'TZ(v-v')^{r',i), (6) 

where the spin-raising kernel 

1 „-iarg(r-r') 



is the Fourier transform of ie~ l ar s( k ). Importantly, the 
7?.-operator in real space has the unitarity property, 
/ drTZ(r — ri)7£(r — r 2 ) = <5(ri — r 2 ), which is the direct 
consequence of the unitarity of the spin-rising transfor- 
mation in the momentum space representation. 

The most notable feature of the dispersion relation (0| 
is that its groundstate is degenerate along the circle in 
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FIG. 1. (Color online) Dispersion relation Q of particles with 
Rashba SO coupling. The degenerate groundstate k = ko — 
mv is shown by full circle. 



the momentum space k = ko — mv. As a result the 
many-body groundstate of ./V non-interacting bosons is 
highly degenerate (not so for fermions, though). Indeed, 
any occupation of the states along the groundstate circle 
k = ko yields exactly the same kinetic energy —Nk^jlm. 
Hereafter we shall measure the energy from that value, 
taking it for the origin of the energy axis. It is therefore 
only the interactions, which may break the degeneracy 
and select a true groundstate. The situation is similar to 
a partially filled Landau level in the context of FQHE. 
There too the kinetic energy is fully degenerate and the 
groundstate is solely determined by the interparticle in- 
teractions. 



B. Bosons with Rashba dispersion 

Let us consider N particles with the s-wave short-range 
interactions of the form 



1 N 

H ^ = ^ £ *k - * + wP^] > 

i,3 



(8) 



where go is the spin-isotropic dimensionless interaction 
constant, while <?2 is the spin-anisotropic interaction con- 
stant. 

One can now evaluate interaction energy of certain 
simple TV-body Bose states. One such state is a Bose- 
Einstcin condensate in a single state belonging to a de- 
generate manifold of the single-particle ground states, i.e. 
a state with momentum k, such that k = fco and 7 = — 1, 



v4 0) = ;QV>k,-(ri,Si) 



N 



(9) 



and k states. We shall call it thus time-reversal sym- 
metry broken (TRSB) state. Note that this state does 
not break the parity P, Eq. ([3]). To see it most clearly 
one may choose momentum k direction to be along the x- 
axis. Calculating the expectation value of the interaction 
energy © over TRSB state, one finds 
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2mV 



go 



(10) 



Since the kinetic energy in the state ((9]) is zero, the in- 
teraction energy (fl"U|) coincides with the total one. 

One may consider now a Bose condensate built on a 
coherent superposition of say two states ki and k2 both 
belonging to the degenerate manifold k — ko: 

N 1 

*fl ) =Il7|^ 1 ,-(r i , S i)+^ k2) -(r i , Si )] . (11) 

i— 1 v 

The corresponding interaction (and thus total) energy is 
found to be 
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go H cos — sin — 

y 2 2 2 2 



, (12) 



where <fi = ar g(ki) — arg(k2) is the angle between the 
two states of the degenerate manifold. It is clear that, 
provided the spin-isotropic interaction is repulsive go > 0, 
the only way the state ([TTj) may be energetically more 
favorable than the state (O is if 52 < 0. In the latter case 
the most favorable choice of the two states corresponds 
to <fr = 7r, i.e. k2 = — ki, with the interaction energy 

£?int = (N 2 /2mV) [go + 92 1 'A- Such a state represents 
a spin density wave (SDW) with a uniform total density 
and the two spin components oscillating harmonically out 
of phase. It is symmetric with respect to both time- 
reversal and parity symmetries, Eqs. @, ©, but breaks 
the rotational symmetry. In either of these states the 
kinetic energy is zero. 

It was conjectured^^ that TRSB ^ ( °\ ©, and the 

SDW state ^\ (HJ), are the many-body groundstates 
of the Rasba interacting bosons for 32 > and < 
correspondingly. It was later suggested^! that the tran- 
sition between TRSB and SDW states is shifted towards 
positive <?2 due to a ceratin admixture of coherently oc- 
cupied BCS-like pairs of k and — k states. We notice 
that for both of these states the chemical potential scales 
linearly with the density n = N/V, 



\.l b = dE- mt /dN oc n . 



(13) 



It is instructive to compare this scaling of the bosonic 
chemical potential with the corresponding fermionic one. 



This wavefunction is obviously symmetric with respect 
to the permutation of any two pairs (rj,Sj) O (rj,Sj). 
Such a state is not symmetric under time-reversal trans- 
formation T, Eq. ([3]), because of unequal population of k 



C. Fermions with Rashba dispersion 

Unlike bosons, the non-interacting Fermi gas exhibits 
the unique ground-state. It is given by the rotationally 
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symmetric Fermi sea with the Fermi surface consisting of 
two concentric circles with the radii kp± = k ± 27m/fco, 
hereafter the small density n <C k 2 , is assumed. The 
corresponding Fermi energy measured from the bottom 
of the spectrum is 

ZTTLKq 

Notice the fact that [ip oc n 2 , which is usually the 
feature of ID Fermi gas. Here it happens because the 
single-particle density of states exhibits ID-like behav- 
ior v(e) oc e~ 1//2 close to the bottom of the Rashba cir- 
cle. The interesting observation is that at small enough 
density n < 5o^'o the chemical potential of interacting 
bosons appears to be bigger than that of the free Fermi 
gas with the same dispersion relation Q . Before making 
conclusions from this observation one needs to consider 
the interaction energy (J8j) of the Fermi sea state. 

Spinlcss (or fully spin-polarized) fermions do not in- 
teract through short-range interactions. In our case the 
particles have spin, which is locked to their orbital mo- 
menta. As a result the symmetric Fermi sea, described 
above, contains all spin directions in x-y plane with equal 
weights. Since two fermions with opposite spins interact 
through the short range interaction ([5]), one expects the 
average interaction energy of the symmetric Fermi sea to 
be of the same order as bosonic one oc g^N 2 /mV and 
thus [If oc n similarly to the bosonic condensates. It was 
recently noticed by Berg, Rudner and Kivelsonii that low 
density Rashba fermions may have paramctrically lower 
groundstatc energy if they form a nematic state. 

To motivate the idea, let us consider two-fcrmion 
state as a Slater determinant built on single par- 
ticle states (JHJ) with momenta ki 2 close to the 
Rashba circle V F = [^,-(1*1, si)V>k 2 ,-(r2, s 2 ) - 
V'ki,-(r2, s 2) - 0k 2 ,-( r ii s i)]/v / 2- The interaction energy 
([SJ of such state is given by Ei nt = sin 2 (0/2)(<7o+52)/2m, 
where 4> = ar g(ki) — arg(k2). Therefore the interac- 
tion energy between fermions tends to zero if arg(ki) — > 
arg(k 2 ), essentially because the spins are aligned in this 
limit and fermions with the same spin can not interact 
through the short-range interaction potential. The BRK 
nematic state takes advantage of this observation. 

To describe such a state qualitatively let us imagine 
that the many-body fermionic state ^Sp is constructed 
as a Slater determinant of states with the angular di- 
rections confined to the angular segment of the momen- 
tum (and spin) space of size O <C 2-7T and momenta close 

to the spin-orbit circle ||kj| — fco| < kp , Fig. [5] The 
corresponding Fermi momentum is found from the con- 
dition 2fc£?' ) fco0 = (27r) 2 n. As a result the correspond- 
ing kinetic energy per particle E kin /N cx [kf^ 2 / m oc 
n 2 /(m& 2 k 2 ). On the other hand, the interaction energy 
per particle is E int /N oc (go + g2)n<d 2 /m, with the factor 
2 originating from sin 2 (</>/2) ~ 2 . One can now min- 
imize the sum of kinetic and interaction energy over O 
to find O oc [n/fcg(<7o + g 2 )] 1//4 , which is indeed small as 



long as 

n < n = k 2 (g a + g 2 ) ■ (15) 

With this O the chemical potential of short-range inter- 
acting fermions with Rashba spin-orbit coupling is found 
to be [Lf oc n 3 / 2 y 'go + g^/mko, which is the result of 
BRK^. In the small density limit (fT5j) the latter is big- 
ger than that of non- interacting Rashba fermions (|14l) . 
but is advantageous over bosonic TRSB and SDW states 
(fl~3l) : (j,p < [ip < hb- Notice also the non-analytic de- 
pendence of [if on the interaction strength, indicating 
the non-perturbative nature of this result. We show in 
Section IIII Dl that the composite fermion nematic state 
may be described quantitatively within Hartree-Fock self- 
consistent mean-field treatment, Fig. [2] 




FIG. 2. (Color online) Hartree-Fock Fermi surfaces of BRK 
nematic states with three different densities. The absolute 
value of the momenta are restricted to the vicinity of the 
Rashba circle k = fco, while their angular directions are re- 
stricted to the angle Q oc [n/kg(go + S^)] 1 ^ 4 - 



D. Can bosons be fermions? 

We have arrived thus to the conclusion that at low 
density (|15|) the chemical potential and groundstate en- 
ergy of fermionic many-body system is smaller than the 
corresponding bosonic one. The similar situation seem- 
ingly happens in a ID system of spinless particles with 
short range interaction potential (go/m)S(xi~Xj). There 
too the mean-field treatment of bosons suggests that 
[is ~ gon. On the the other hand, spinless fermions 
(which are not affected by short-range interactions due to 
the Pauli principle) exhibit [ip ^ n 2 . It thus seems that 
at a small enough density the fermionic groundstate en- 
ergy is smaller than the bosonic one. The actual situation 
is, of course, very different. At small density bosonic 
many-body groundstate wavefunction ^ b approaches the 
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FIG. 3. (Color online) Equations of state /x(n) of ID quantum 
gases with short-range interactions: spinless Fermi gas (long- 
dashed); mean-field approximation for Bose gas (dashed); ex- 
act result^ for bosons (full). Notice that bosonic ground- 
state energy E = / /J,dn is always smaller than the fermionic 
one, although the mean-field treatment suggests otherwise for 
n < no = go- 

symmetrized fermionic one 

N 

&b(xi, ■ ■ -,Xn) = JJsign^i - Xj)^f F (xi,. . .,x N ) , 

(16) 

where $f is the fermionic Slater determinant occupying 
a finite portion of the momentum space — im < k < im. 
It is important to notice that, although sign(xi — Xj) is 
undefined at x% = Xj, the fermionic part cancels at all 
such points making fg^!,...,!^) well-defined in the 
entire space. Due to the same observation there is no in- 
teraction energy cost for the short-range repulsion. The 
corresponding kinetic energy per particle oc (irn) 2 /2m is 
small in the low density limit. This is the so-called Tonks- 
Girardeau limi t 38 ' 39 i , where bosons redress themselves 
as fermions. This allows them to take advantage of the 
wavefunction, which is nullified at all points where any 
two particles approach each other. As a result they avoid 
paying short-range interaction energy cost, while corre- 
sponding kinetic energy cost appears to be worth the 
bargain at low density. The corresponding ID equations 
of state are depicted in Fig. [3] 

It is thus tempting to speculate that our spinfull 2D 
system may benefit from a similar construction. Namely, 
at low density bosons may want to redress themselves 
as fermions to take advantage of their lower interac- 
tion energy. The blueprints of the boson-fcrmion cor- 
respondence in 2D are provided by FQHE studie o 45 ' 47 , 
where the correspondence is achieved by ascribing Chern- 
Simons phase factor to many-body wavefunctions. Such 
factor takes care of the proper symmetry of wavefunc- 
tions. It comes with the price however: the gauge mag- 
netic field which has a form of delta-functional flux tubes 
attached to every particle. It is believed that (under 
proper conditions) this magnetic field may be substi- 



tuted by a uniform one, making the problem analytically 
tractable. In the next section we develop a similar strat- 
egy for the chiral spinfull boson-fcrmion correspondence. 

III. COMPOSITE FERMION STATE 

A. Naive fermionization attempt 

Our goal is to construct a variational groundstate for 
Rashba bosons based on the composite fermion idea. We 
argue here that in the small density limit such a state, 
inspired by the Tonks-Girardcau limit (|16[) . is advanta- 
geous over both TRSB © and SDW (HI]) bosonic states. 
The two main differences with the Tonks-Girardcau case 
are: (a) the 2D nature of our problem and (b) the spinfull 
nature of the particles. It is still tempting to straight- 
forwardly generalize Eq. p^|) to the 2D spinfull case by 
taking the variational ground state of the Bose gas in the 
following form: 

* B (ri,si,---rjv,Sjv) (17) 
= Y[e iaT ^-^* F (r 1 ,s 1 ,---T N ,s N ), 

where ^/^(ri, s\, ■ ■ ■ r^f , sjv) is an N-particle fermionic 
wavefunction. The Chern-Simons phase e 4arg ( ri ~ rj ) = 
(zi — Zj)/\zi — Zj\, with Zj labeling complex spatial co- 
ordinates, Zj = Xj + iyj, is antisymmetric with respect 
to exchange of any two coordinates. As a result, many- 
body wavefunction, ^s(ri, si, ■ ■ ■ r^r, sat), is symmetric 
under permutation of pairs , Si and Tj , sj of coordi- 
nates and spins, Sj =j" , ^, if the fermionic wave function, 
^^(ri, si, • • • r^v, sn), is antisymmetric with respect to 
the same permutations. One way of writing the latter is 
to take it as N x N Slater determinant of e.g. single- 
particle spinor wave functions, ^kj,-(i"i, Si), Eq. ([5]), 
where i, j — 1, 2, . . . , N. 

While being probably the most straightforward way of 
guessing a spinfull composite fermion state, Eq. (fTT|) has 
a number of fatal drawbacks. First, one might naively 
expect that similarly to Eq. (|16[) , this ansatz maps the in- 
teracting bosonic system onto a system of non- interacting 
fermions. A closer look shows that this is not the case. 
Indeed, although the Slater determinant implies that the 
fermionic wave function has zeros at coinciding spatial 
points and spins, however for coinciding points and dif- 
ferent spins the wavefunction is not nullified. As a result 
the composite fermions with opposite spins still do in- 
teract, despite the short-range nature of the interaction 
potential. 

Even more serious problem with the trial wavefunc- 
tion (JTTJ) is that it is not well-defined for coinciding 
points and opposite spins. This is due to the fact that 
the Chern-Simons phase is singular at coinciding points, 
while the fermionic part typ is not vanishing if spins are 
opposite. This ambiguity leads to logarithmic divergent 
contributions to the average kinetic energy of the state 
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(fT7f . Indeed, consider the part of the kinetic energy 
(2m) _1 J drilVr^sl 2 , where the gradient operators act 
on the Chern- Simons phases. This leads to the kinetic 
energy contribution of the form 



1 

2m 



|* F (ri,si, • ■ -rjv, s N ) 



(18) 

The diagonal terms j = j' in this double sum lead to the 
integrals of the form J dri/|rj — r j | 2 , which exhibit loga- 
rithmic behavior when r< w r 3 . If particles i and j have 
the same spin, typ = at r, = Fj cutting the logarithmic 
divergence of the integral at small distances ~ |kj — kj| _1 
(at large distances it is cut at a typical interparticle dis- 
tance n -1 / 2 due to the random sign of the numerator in 
Eq. (|T8|) ). However for opposite spins &p ^ at = r_, 
and the integral (jlijp diverges in all such points. The re- 
sult is logarithmical divergent chemical potential, mak- 
ing the trial wavefunction (|T71) essentially useless. One 
should thus look for an alternative way to introduce com- 
posite fermion state for spinfull Rashba particles, which 
avoids logarithmic divergent terms in the kinetic energy. 



B. Fermionization 

The idea for an alternative scheme comes from the ob- 
servation that at small density all relevant energy scales 
are much smaller than SO energy eo = mv 2 /2. There- 
fore one would like to have a many-body state which is 
projected onto the Hilbcrt subspacc of the lower spin- 
orbit branch 7 = —1, Eq. (£[]). On the single particle 
level such a projection is achieved by ensuring the rela- 
tion between up and down components of the spinor. 
One can straightforwardly generalize it for a many-body 
wavefunction. To this end one should specify a fully sym- 
metric in the coordinate space wavefunction of the min- 
imal spin 



*4..„4.(ri, ...,T N ) = *(ri |, . . . , r N 1) . 



(19) 



Then all other spin components of the fully projected 
wavefunction may be uniquely determined from the min- 
imal spin component by successive application of the spin 
raising non-local operator 1Z, Eq. ([7]), e.g. 



*(nt, 



,r N i) =jdv' l n{T 1 - ri)* 4 .. 4 (ri, . . . ,r N ) , 



tt(n t, . . . , rjv t) = / dr[ . . . dr' N K{v 1 - r'j) . . . K(r N - r' N ) 



(20) 



It is easy to see that all the components defined this way 
are symmetric with respect to simultaneous interchange 
of r j , Si and r j , Sj . 

One can now fermionize such a wavefunction by writ- 
ing the spatially symmetric minimal spin component (|19[) 
as a product of Chern-Simons phase and fully antisym- 
metric spinless fermionic wavefunction. The latter will 



be shown to describe the anisotropic nematic state. It is 
therefore convenient to incorporate the same anisotropy 
in the Chern-Simons phase too. We thus define rescaled 
coordinates x = ax, y = y/a and f = (x, y), where the 
density-dependent scaling parameter a will be specified 
in Section IIII El The fcrmionized minimal spin compo- 
nent is then written as: 



*|...i(ri, 



1 



)JV 



(21) 

where we have introduced chirality factor A = ±1, which 
defines the direction of the Chern-Simons flux relative 
to the spin chirality. The many-body fermionic wave- 
function f(ji, . . . , rjv) is antisymmetric with respect to 
permutation of any two of its spatial arguments. Thanks 
to the Chern-Simons phase factor the left hand side of 
Eq. (|2ip is fully symmetric both in coordinate and spin 
spaces. Notice that, unlike the earlier attempt, Eq. (|17l) . 
the wavefunction ([2"T]) is everywhere well-defined. This 
is the case because the Chern-Simons phase multiplies 
the fully antisymmetric function of spinless fcrmions, 
which cancels if any of its spatial arguments coincide. 
This became possible by adding the Chern-Simons fac- 
tor to a fully spin-symmetric minimal spin component 
only. The higher spin components are then built up from 
the fermionized minimal spin bosonic function (|21[) by a 
successive application of the spin raising kernel, Eq. (|20[) . 
As a result all 2 N spin components of the TV-body bosonic 
wavefunction (|2T)|) . (f2~Tj) are well-defined functions in the 
entire 2iV-dimensional coordinate space. This is unlike 
the earlier attempt, see Eq. (15), which resulted in the 
logarithmic divergences in the kinetic energy. 

Of course, there is a complimentary construction 
where one fermionizes the maximal spin component 
^t"'t( r i> • ■ • ' Tn ) an d builds the lower spin components 
by successive action of spin lowering operator. These 
two states are related by the parity transformation ([3]) 
acting on the coordinates and spins of all N particles. 
Namely, acting with such parity operator on the state 
(|20)) . (f2~Tj) built from the minimal spin with chirality A, 
one obtains the state built from the maximal spin compo- 
nent with chirality — A and parity transformed, yj — > —yj, 
fermionic wavefunction Indeed, the parity operation 



Q interchanges spins and transforms Zj 



resulting 



in A — > —A correspondence. Since the parity operator P 
commutes with the full Hamiltonian Hq + H- mt , these two 
states are degenerate. 

Let us now discuss the average kinetic energy of the 
state with the wavefunction (flU|) . (|2TH) . Evaluating the 
expectation value of the single-particle operator incorpo- 
rating kinetic and spin-orbit parts ([1]), one finds 

„N N 

Emu = 2 N m dr^ ^(n, . . . , v N ) £ K(r 3 - rj) 

x^.. 4 (r' 1 ,...,r' JV ), (22) 
where the non-local operator of the kinetic energy acting 
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on the minimal spin component is 



K(r-r') = — 



2m 



--[^(r'-r^-^r-r')], 

(23) 

here we employed notations df- = d x ± id y . This form 
is obtained by expressing the spin-up states in terms of 
the spin-down states with the help of Eq. ([6]) and us- 
ing unitarity of the spin-raising operator. The Fourier 
transformation of the kernel (|23|) results in the lower 
branch, 7 = — 1, of the single particle spectrum Eq. (|4]) 
K{k) = k 2 /2m — u|k| = which is, of course, ex- 

pected for the projected wavefunction in the form (|19l) . 
([20|) . It is exactly the reason to build the higher spin 
states with the help of the spin-raising operator ([2"U| to 
ensure that the kinetic energy belongs to the lower spin- 
orbit branch. 

The non-analytic behavior of the spectrum (0| at k = 
translates into the non-local behavior of the kinetic en- 
ergy kernel (|23|l in the coordinate representation. This 
non-locality complicates the way the kinetic energy oper- 
ator acts on the Chern-Simons phase in Eq. (|2"Tj) . On the 
other hand, the low energy part of the Hilbert space is 
located close to the Rashba circle |k| = fco, i.e. far away 
from the k — singularity. Below we discuss variational 
choices for the many-body fcrmionic wavefunction typ, 
which explicitly include only momentum components lo- 
calized around |k| = fco circle. For those components 
the non-analyticity at fc = and thus non-locality of 
the kinetic energy in the coordinate space are not essen- 
tial. Therefore the single-particle kinetic energy spec- 
trum near its minimum may be approximated by an an- 
alytic function of momentum (and thus local differential 
operator in the coordinate representation) as, e.g., 



K(k) = -co 



(|k| 2 -fc 2 ) 2 /(8mfc 2 ). 



(24) 



Its action on Chern-Simons phase results into the substi- 
tution of the gradient operators by the covariant deriva- 
tives with the gauge vector potential (see below). 

The important observation at this stage is that there 
are no logarithmical divergent contributions to the ki- 
netic energy, which were fatal for the earlier fcrmioniza- 
tion attempt (|T7|) . Indeed, the kinetic energy kernel (|22|) 
acting on the Chern-Simons phase ([2"Tj) results, among 
others, in the term similar to (fT5]l . However, this time the 
fermionic wavefunction p does not contain spin indices 
and cancels any time t% = Tj. As a result, all the loga- 
rithmic integrals are cut off at small distances by a scale 
built into the fcrmionic wavefunction ^ p(ri, . . . , rjv)- At 
large distances they are cut off at a typical intcrparticle 
distance. As we discuss in Section UlI El these two length 
scales can be made of the same order by an appropriate 
choice of the rescaling parameter a in the Chern-Simons 
phase (f2~Tj) . This makes the logarithm to be a number of 
order one, showing that the localized flux-tube structure 
of the Chern-Simons magnetic field (see below) does not 
strongly affect the kinetic energy of the fermion state fyp . 
It suggests the mean-field substitution of the flux tubes 



by a uniform magnetic field, analogous to those employed 
in FQHE literature, e.gM^L. 

Next we discuss inter-particle interactions, excluding 
for a while the Chern-Simons term in the wavefunction 
(f2"l"j) . Following BRK, we employ Hartree-Fock approxi- 
mation to minimize the total energy of interacting par- 
ticles having spectrum (|2"4")) and find that the chemical 
potential scales as [i ~ r?l 2 . Finally, we include the 
Chern-Simons term within the Hartree-Fock approxima- 
tion for interactions. We argue that it does not affect 
the parametric dependence oc n 3 / 2 , although does 
affect the spectrum of single-particle excitations. 



C. Interactions 

We focus now on the average energy of the short- 
range interactions ([5]) in the many-body state specified 
by Eqs. (|2"Tj) and ([2H)l . Since the minimal spin compo- 
nent (|2T]) cancels at coinciding spatial points, it does not 
contribute to the interaction energy. On the other hand, 
the higher spin components, built according to Eq. (|2"0|) . 
do not vanish at coinciding points and thus lead to a 
non-zero interaction energy. At the first glance this ob- 
servation makes the interaction energy of the composite 
fermion state Eqs. (|2"Tj) and (j2T)]) to be of the same or- 
der as that of the bosonic condensate (|9|) or (fTTj) . mak- 
ing the entire construction useless. We show below that 
this is not the case thanks to the antisymmetric nature 
of ^/p(ri, . . . ,rjy). The latter leads to the Hartree mi- 
nus Fock structure of the two-particle interactions, which 
nearly cancels the interaction energy for particles with 
nearly collincar momenta (and thus spin) directions. On 
the other hand, the bosonic condensate © admits only 
the Hartree contribution, while ([TT]) leads to Hartree plus 
Fock structure, not exhibiting the cancelation. 

To derive the effective two-body interactions in the 
composite fermion state it is sufficient to consider two 
particles in such state built of two single-particle states, 
e.g. ^kj(r) = e 4kj " r /\/y with j = 1,2 as 



*Hri,r 2 ) 
*u(ri,r 2 ) 



1 

1 



■yMi-i^ksfo) - ^(^i/Mri) 



- _ p *Aarg(ri— Pa) 

~ 2 



*F(ri,r 2 ), 



(25) 



*n(ri,r 2 ) = J drift(n-ri)¥ u (ri,r 2 ), 
* it (n,r 2 ) = J dr 2 7e(r 2 -r 2 )* ;i (n,r 2 ), 

* n (n,r 2 ) = J dv^n^-v^n^-v'^^v'^). 

The average interaction energy according to Eq. ([S| is 
given by 



L 1 



(. 9o+52 )|vI/| A t )(r,r)| 



(26) 



+ (ffo - ff 3 )|*g ) (r,r)| a + ( 9o fl2 )|^(r,r)| 
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and therefore one needs to evaluate the wavefunction 
for higher spin components at coinciding spatial points. 
They are given by: 



ir(ki+k 2 ) 



Ft 



ki,k 2 



*Tt( r > r ) = ^7 e 
* n (r, r) = *^(r, r) = -L e *<*+><») Gkl>1 



(27) 



where ^ki, k 2 an( i G , k ll k 2 are interaction form- factors eval- 
uated in the Appendix. In the approximation where 
|ki| ps [lea [ ps fco and arg(k^ 2 ) <C 27r we found: 



■Fk!,k 2 = ic(arg(ki) - arg(k 2 )) , 
G kl ,k 2 = d(axg(ki) - arg(k 2 )) . 



(28) 



Here c = c\(a) and d = d\(a) are numerical factors 
weakly dependent on CS chirality A and anisotropy a. 
For example, |c+(l)| ~ 1.22; |c_(l)| ~ 0.85; c±(0) = 
and d+(l) = 0; ^ 1.41; |d±(0)| ~ 0.90. A 

very important observation is that the interaction en- 
ergy tends to zero if arg(ki) — > arg(k 2 ). The relative 
minus sign in Eqs. (|28p is entirely due to the compos- 
ite fermion nature of the wavefunction (J25J). Indeed 
the two-particle fermionic wavefunction (|25[) is zero if 
arg(ki) = arg(k 2 ) due to Pauli blocking. This offers 
a possibility to construct a trial many-body wavefunc- 
tion which takes advantage of the fact that the particles 
with nearly collinear spins (and thus momenta!) inter- 
act only weakly. As discussed in section III C[ the similar 
construction for fermions with the Rashba SO coupling 
was recently employed by BRK— , who showed that the 
optimal trial state is a nematic one. 



D. Hartree-Fock theory 

The secondary quantized interaction Hamiltonian for 
projected spinless fermions takes the form 



1 



E 



M, 



ki,k 2 J J 



32m ^ " ik i' k = c k 2 c ki c kiCk 2 

k' 1 +k2=ki+k 2 



(29) 



where the interaction matrix elements 



(go + 52 )i r k' 1 -Pki ,k 2 + 2 (50 — 52 )Gk[ ,ki, ,k 2 



M k /' k , 2 

(30) 

depend on all four (two incoming and two outgoing) mo- 
menta and therefore interactions can't be reduced to the 
density-density form. Moreover in the coordinate repre- 
sentation such interactions acquire essentially non-local 
|r — r'|~ 2 form, which is a consequence of the projec- 
tion on the lower SO branch^. In principle the projec- 
tion, Eqs. (|19p and (|20[) . generates also three- and more- 
particle interactions. One may check however that in the 
dilute limit ()15[) their effect is negligibly small in com- 
parison with the two-particle term kept in Eq. (|29[) . 

One can now treat the interaction term in Eq. (|29|) 
in the Hartree-Fock approximation by pairing one cre- 
ation and one annihilation fermionic operator (c k ,Ck) = 



(a) k, 



o 



k (b) k 



k (c) 




FIG. 4. (Color online) (a) Two-body composite fermion inter- 
action with the matrix element M. /'. ?, Eq. (1301) . (b) Hartree 

and (c) Fock ways of pairing the operators in the interaction 
Hamiltonian (1291). 



<5(k' — k)rik'. The Hartree and Fock ways of pairing are 
illustrated in Fig. H^,b, correspondingly. This way one 
obtains the mean-held single-particle Hamiltonian 



Hn F = ^2 i^HF(k) c k c k , 



where the Hartree-Fock kinetic energy is given by 



(31) 



k.k' 
k,k' 



k',k 
k.k' 



n k ', (32) 



and the self-consistency condition is imposed by requiring 
?l k' = /((A'(k') — hf)/T) is the equilibrium occupation 
number of the state k' determined by its total energy 
A'(k') and the chemical potential \ip. The latter is to be 
found from particle conservation J^ k , = n. 

In the limit of small density (|15[) we anticipate the 
nematic state, discussed qualitatively in section III CI 
Choosing its center to be along the positive a;-direction, 
one may approximate arg(k) ps k y /ko -C 27r. We can now 
employ the form-factors ([28]) as well as the lower branch 
dispersion £k,_ in the form of Eq. (|2~4"|) . to hnd for the 
Hartrcc-Fock energy 



A HF (k) « E Q 



(k x - fc ) 2 
2m 



1.2 

K y 
2m, 



where 



Akl 

m v = m ^> m 

gn 



(33) 



(34) 



and E = g ^2k'(ky) 2n k' / (4mfcp) with the effective inter- 
action constant g = (go + .92) |c| 2 + 2(go — 92)\d\ 2 . We 
note here that m y is essentially insensitive to the distri- 
bution rife, as long as it covers a small fraction of the 
Rashba circle in the momentum space. It is clear that 
the Fermi surface, corresponding to the dispersion rela- 
tion (f3"3"l) is an ellipse elongated along the y direction and 
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centered around (fcojO), Fig. [2j One can now find the 
chemical potential of interacting fermions at T = by 
solving J2k(]i)<h f ~ n - This way we find for the chemi- 
cal potential of the interacting composite fermions 



5ir n 3 / 2 



(35) 



and Eq = /Uf/5, in agreement with the estimate in the 
end of Section IIICl 



E. Chern-Simons gauge field 

So far we have been discussing the interaction energy of 
the composite fermions. The latter are related to bosons 
through the Chern-Simons phase according to Eq. ([21]) . 
We need to discuss now the role of this phase. At the 
Hartree-Fock level the system is effectively described by 
non-interacting quasiparticles with the anisotropic dis- 
persion relation Eq. ()33|) . In the coordinate representa- 
tion the momentum operators are given by k x , y — > id x , v - 
To bring it into a more familiar form one may shift the 
momentum origin to (ko, 0) by the gauge factor e tk ° x and 
rescale the variables as x = ax and y = y/ot, where 
a = (m/rriy) 1 / 4 . In the rescaled coordinates the Hartree- 
Fock Hamiltonian is isotropic with the cyclotron mass 
m c = ^jmniy. Upon acting on the Chern-Simons phase 
the rescaled momentum operator results in a vector po- 
tential A(f), which is included by replacing components 
of iVf by the corresponding components of iV? — A(f), 



^HF = J2 i - A *(^)) 2 + ( id V> - A ^)) 2 ] , 

(36) 



3=1 



where the vector potential is given by 5 -^ 

( f 3 - f 0/3 



A «( f j) = J2 e 



a/3 



W-i - Vi 



while the corresponding Chern-Simons magnetic field di- 
rected perpendicular to the 2D plane is given by 

Bfa) = curlAfa) = 2tt^ Sfa -f*) = 2tt^ tffo -r t ). 

(38) 

To obtain the isotropic form of the Hamiltonian acting 
in the fermionic space, it is important that the Chern- 
Simons phase was defined with the rescaled coordinates 
f, Eq. (|2ip . We can now specify the value of the rescal- 
ing parameter a = (m/my) 1 / 4 cx (n/?io) 1//4 < 1. Notice 
that since m y is itself weakly dependent on a (through 
interaction form-factors c an d) the above definition of 
a is actually a self-consistent equation. Moreover, in- 
clusion of the magnetic field affects the wavefunctions of 
the quasiparticles, e.g. a homogeneous field results in the 
Landau quantization. That, in turn, modifies the matrix 
elements of interactions M^}'^? defining m y . However, 



because of the form of interaction -Fki,k 2 and Gki,k 2 the 
effective mass appears to depend only weakly on the spe- 
cific form of the wave functions, as long as those com- 
posed of plane waves with wave vectors in the vicinity of 
the point (fco,0). This appears to be the case even for a 
quantizing magnetic field corresponding to a fully occu- 
pied single Landau level. Therefore we use Eq. (|34| for 
m y in the following and apply the effective description 
Eq. (|36p even for a quantizing magnetic field. 

The importance of bringing the effective fermionic 
Hamiltonian to the isotropic form (|36|) is to argue that 
the kinetic energy does not contain large logarithms. As 
explained below Eq. fjl8j> . the kinetic energy in presence 
of the singular magnetic field (|38p contain logarithmic 
integrals. At small distance they are cut by the scale of 
the correlation hole in the fermionic wavefunction Vfj?, 
while at large distance they are cut by the average dis- 
tance between the particles. In the isotropic fermionic 
state described by the Hamiltonian (|36[) both of these 



two length scales are given by kp 1 ~ n^ 1 / 2 . As a result, 
the logarithms contain a number rather than a density- 
dependent parameter. Notice that without rescaling of 
CS phase such logarithms would contain log(l/a), mak- 
ing the kinetic energy of the order n 3 ! 2 log(no / n) . Rescal- 
ing avoids this large extra factor in the kinetic energy. 

Keeping the kinetic energy to be cx r?l 2 by the 
fermionic nature of the state along with the careful 
rescaling of CS phase, suggests to employ the mean-field 
treatmen t 45 ' 47 of CS magnetic field. It substitutes the 
collection of flux lines (|38p by a uniform magnetic field 
with the same total flux. The latter is given by 2D den- 
sity of particles (it is important that the rescaling of CS 
phase is area-preserving and thus not affecting the total 
flux), 



B = 2im . 



(39) 



The corresponding mean-field vector potential may be 



(37) chosen as A x (rj) = —nnijj and A y (fj) = 



: P In 

this approximation the Hamiltonian (|36p represents a 
one-body problem of a particle with the cyclotron mass 
m c = ^/mrriy cx m-^rio/n in a constant magnetic field 
(|39p . The corresponding Landau spectrum is e; = 
u> c (l + 1/2), where I = 0,1,... and the cyclotron fre- 
quency uj c = 2im/m c . Since by construction we intro- 
duced exactly one flux quantum per particle, the corre- 
sponding filling factor is v = 1. The A-body ground- 
state wavefunction ^^(f i, . . . , fjy) is thus given by the 
Slater determinant built from the single-particle wave- 
functions of the fully occupied lowest Landau level (LLL) 
(the center of mass shift in the direction of the ne- 
matic order produces an additional multiplicative factor 
exp{iko 53 - Xj}). In terms of rescaled complex coordi- 
nates Zj = Xj + ijjj = axj + iyj/a the corresponding 
fermionic wavefunction takes the fori 



,50. 



A' 



i<3 



-irn J2 \Zj | 2 /2 ikoJ2 x 3 



(40) 
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where Cat is normalization factor. The Chcrn-Simons 
phase may be also written in terms of the rescaled com- 
plex coordinates as Y\{zi — Zj)/\zi ~ Zj\. As a result 
the minimal-spin component of the bosonic wavefunction 
pTj) takes the simple form 



N 



- N,2 c N \{ 



-7T77 Yl 



Zi 



2 /2 ik 

e 



(41) 

The higher spin components are obtained by acting with 
the spin raising operators on the minimal spin compo- 
nent, Eq. (f2"D")) . Notice that the spin-raising operator ([7]) 
is to be written in non-rescaled original coordinates Vj to 
ensure that the wavefunction is projected on the lower 
SO branch. The wavefunction (|4Tj) is independent of 
the CS chirality A. It depends on the single density- 



1/4 



oc 



and interaction-dependent parameter a — (m/m y ) 
(ji/jio) 1 / 4 , which specifies the anisotropy. This wavefunc- 
tion describes an ellipsoidal droplet of the gas elongated 
along the x direction with the ratio of x and y axes given 
by a 2 . The average density in this droplet is n, while its 
size depends on the number of particles N. The state 
(j4Tj) clearly breaks rotation symmetry. It also breaks 
time-reversal symmetry as well as parity. Indeed, the 
parity operation transforms it into the state descendant 
from the maximal spin component, which is clearly a dif- 
ferent, degenerate state. The average energy per particle 
for the Hamiltonian (|3l))) in such variational A-body state 
is given by 



fip — E + e 



3tt n 3 / 2 
4 mho 



(42) 



It indeed confirms the expectation that the composite 
fcrmion function of Eqs. (|2T))) and (|21[) yields the energy 
lower than TRSB and SDW bosonic states. 



IV. DISCUSSION OF THE RESULTS 

A. Phase Diagram 

As we have seen, the composite fermion energy per 
particle is /if ~ HBy/n/n-o, where (j,b, Eq. (fT3|) . is the 
chemical potential of bosonic TRSB and SDW states. 



Therefore at small density n 



< 



gk^ the compos- 



ite fermion (CF) groundstate is energetically favorable 
over bosonic states. On the other hand, at larger density 
n ^ n o the bosonic states, discussed earlie r 15 ' 23 " —, have 
lower energy. We expect thus to observe quantum phase 
transitions as functions of the chemical potential \i as well 
as anisotropy of the interactions g^. The corresponding 
phase diagram is schematically depicted in Fig. [5l 

To find out about the nature of the transitions we recall 
that all three phases break rotational symmetry. In ad- 
dition TRSB breaks time-reversal, while CF breaks both 
time-reversal and parity symmetries. We expect thus the 
first order transition between CF and SDW, where the 



TRSB 




FIG. 5. (Color online) Phase diagram of Rashba SO bosons. 
The full lines are first order transitions. Insets: schematic 
representation of occupation factors in momentum space (red 
dots are coherent bosonic condensates, blue ellipse is CF 
Fermi sea). 



two discrete symmetries t and P, Eqs. ©, 0, are bro- 
ken simultaneously. Also the transition between SDW 
and TRSB, taking place close to 52 = 0, is to be of the 
first order. Indeed, to avoid density wave modulation in 
SDW phase the populations of two opposite states on the 
Rashba circle must be exactly equal. Therefore the tran- 
sition into TRSB phase with only single populated state 
must be of the first order. 

Transition from TRSB into CF phase breaks the par- 
ity symmetry and could be of the first or second order. 
One can investigate stability of the bosonic TRSB state 
against introduction of the small CF component. To this 
end one can write a variational wavefunction, which con- 
tains superposition of TRSB condensate, Eq. with 
Nb particles and nematic CF fraction with Nf <C Nf, 
particles and Nb + Nf = N, 

*(ri, . . . ,rjv) cc^2 *s ( r i> • • ^ r N b )^CF{r Nb+ i, . . .,r N ), 

(43) 

where P stays for a permutation of arguments between 
Nf, and Nf sets and we have suppressed spin indices for 
brevity. It is important to realize that the bosonic con- 
densate and a small nematic CF fraction prefer to be on 
the opposite sides of the Rashba circle, Fig.[SJ Indeed this 
way the spin parts of respective wavefunctions are (al- 
most) orthogonal, minimizing the exchange energy. The 
residual exchange interactions, due incomplete orthogo- 
nality, lead to CF effective mass m y in the direction tan- 
gential to the Rashba circle. The energy per unit volume 
for the state (|43[) is given by 



E = 



Sirnj („ - n/) 2 



9o 



(n - n f )n } 



[.90 + 52], (44) 



4m c 2m m 

where n/ = Nf/N and the fermionic cyclotron mass 
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FIG. 6. (Color online) (a) Specific energy vs. density for 
CF and Bose condensate phases. The two graphs intersects 
at the critical density n c . In the range hf < n < ub there 
is the phase separation, (b) Chemical potential vs. density. 
Maxwell construction is shown. Compare it to ID case, Fig. [3] 



m c = ^Jmm y = mko / y/ (go — gi)nj2 originates from 
the exchange interactions of the CF fraction with the 
bosonic condensate (in the limit rif <^ n the interactions 
between composite fermions are less important). This 
mean-field expression shows that for gi > TRSB state, 
i.e. rif = 0, is always the energy minimum and thus it is 
stable against small CF fraction. An additional energy 
minimum develops at n/ = n for n < 2no/37r. Although 
one should not take Eq. (|44|) as quantitatively accurate 
beyond n/ <C n limit, it is true that at small enough 
density the CF phase costs less energy than the bosonic 
condensates. Together with the local stability of the con- 
densate it suggests the first order transition into the CF 
phase. 

For the first order transition the equation of state im- 
plies the range of density where the Bose condensate 
(TRSW or SDW) fraction separates from CF fraction. 
This is illustrated on Fig. [6^- In the region rip < n < Ub, 
determined by the common tangential to Ep(n) and 



Eb{u), it is energetically favorable to spatially separate 
the Bose condensate with the fixed density from the 
CF liquid with fixed density rip. Changing the total den- 
sity within this window results in the change between the 
relative fraction of the two, not changing their individual 
densities. In terms of the chemical potential vs. den- 
sity it means a constant critical chemical potential /i c , 
found from the Maxwell construction, within the window 
np < n < riB, Fig.[6}D. Alternatively, it means a range of 
a constant critical quantum pressure p c , Fig. |6K, on the 
pressure vs. volume T — isotherm. To minimize the 
interfacial energy TRSB and CF fractions prefer to have 
opposite momenta and spin. This allows avoiding paying 
the exchange interaction energy thanks to orthogonality 
of the spinors. One expects thus to find antiferromag- 
netically ordered mixture of the TRSB Bose condensate 
and CF fractions. 



B. Composite fermion state in a trap 

Consider an axially symmetric trap created by a har- 
monic potential V(r) = ttuJqT 1 . In the Thomas- Fermi 
approximation the local density n(r) is found from the 
condition V(r) + fx(n(r)) — fi = 0, where fx(n) repre- 
sents microscopic equation of state and fi is the macro- 
scopic chemical potential, found from the condition 

2ir J^ 1 rdrn(r) = N, where V(R) — \i. For CF equation 
of state \x(n) oc ^fgr?!" 1 /mko, this program yields: 

(45) 

Here R represents the spatial extent of the Af-particle 
groundstate, while /i is the typical kinetic energy per par- 
ticle measured upon trap release. These results should 
be compared with those for the bosonic condensate with 
the equation of state /u(n) cx gn. The latter yields^ 
/i cx uj^gN) 1 / 2 and R oc N 1 / 4 . The CF scaling of the 
chemical potential in the harmonic trap fip oc N 3 ^ 5 is 
valid as long as it less than the corresponding bosonic 
result ub oc A 1 / 2 . Equating the two of them, one finds 
the condition for the number of particles within the trap 
below which the CF groundstate prevails 

N<No=g "{^w) 2 ■ (46) 

This is exactly the condition of having the density in the 
middle of the trap less than the critical one n(0) < n c . 
The density profile acquires the shape n{r) oc (i? 2 — 
r 2 ) 2 / 3 . As opposed to the Bose condensate profile^ 3 - 
n(r) oc (R 2 — r 2 ), it exhibits infinite slop at the outer 
edge r ps R. 

In experiments of Refs. [l[ and [|| synthetic gauge field 
and SO coupling had been engineered with the help of 
A = 804. lnm Raman lasers, leading to the typical SO 
momentum fco = 2ir/\ and energy scale k^/m ~ 7kHz. 
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Taking a trap frequency ujq — 30Hz as in e.g. Ref. [551 ] . 
one obtains fc 2 /mw ~ 0.23 x 10 3 and N ~ g 3 x 10 5 . 
The effective coupling constant g in 2D gas of Rb with 
a s = 55A was estimated to be 5 -^ g ~ 0.2, leading to 
Nq ~ 10 3 , which is achievable with modern detection 
techniques. Interestingly enough, studying much smaller 
particle number experimentally became feasible with the 
development of single-atom detection technolog y 5 ^ 58 . 

Since the CF phase has filling factor one, we expect it 
to be gaped in the bulk of the trap. At the surface it 
supports a chiral edge mode, which is similar to v = 1 
quantum Hall edge. In this sense CF state of Rashba SO 
bosons is an interacting topological insulator. 

At a larger particle number N > No the middle of the 
trap has the density exceeding the critical one and thus 
it reverts to one of the Bosc states (TRSB or SDW). The 
density decays towards the edges of the trap, reaching ns 
at some radius. Here the phase separation, Fig. |H1 takes 
place and the density discontinuously drops down to uf ■ 
The outer periphery of the trap contains "vaporized" low- 
density CF phase, while its inner core is filled with the 
"liquid" high-density Bose condensate phase. The overall 
scaling of the chemical potential with the particle number 
approaches the Bose one, fj, ~ iV 1 / 2 . 



C. Rotating systems 

The time-reversal and parity broken CF state is chi- 
ral. It is not immediately obvious from the minimal 
spin component wavefunction (|4Tj) . since it depends on 
the absolute values |5, — Zj\ only. However, the spin- 
raising kernel (J7]) is certainly chiral and so are all the 
higher spin components of the many-body bosonic wave- 
function (|20[) . The degenerate state descending from the 
maximal spin component (which has exactly the same 
form as Eq. (f4"Tj) ) has the opposite chirality. The CF 
groundstate spontaneously breaks the Zi symmetry be- 
tween them. Rotation of the system serves as an ex- 
plicit symmetry breaking perturbation, enforcing one of 
the sates vs. another. Indeed, rotation with the angu- 
lar frequency Q may be viewed as an external magnetic 
field 5 ^ B rot = 2m c f2z, which adds up to CS magnetic 
held B = ±27mz. As a result the total flux is bigger or 
smaller than one flux quantum per particle, depending 
on whether CS magnetic filed is parallel or antiparallel 
with the rotation direction. This either creates m c f2/7r 
holes per unit area in LLL, or promotes the same num- 
ber of particles to the next Landau level. At least from 
the standpoint of the kinetic energy the former alterna- 
tive requires twice less energy than the latter. One ex- 
pects thus that the symmetry is broken in a way to add 
constructively CS and rotational fields to create holes in 
LLL. Due to presence of these holes, one expects gapless 
bulk excitations in the rotating system. (Alternatively in 
analogy with FQHE, excitations may have gaps substan- 
tially reduced compared to the ft = case.) 



D. Fractional Hall states? 

We have employed Chcrn-Simons phase with one flux 
quantum per particle to convert composite fcrmions into 
bosons, i.e. A = ±1 in Eq. (|21l) . One could also attach 
three flux quanta by choosing A = ±3. This choice leads 
to the effective magnetic field B = Girn and thus to v = 
1/3 filling factor. The CF variational groundstate is then 
given by the Laughlin 5 ^ state, i.e. (£; — Zj) — > (£; — zj) 3 
in Eq. (|^0|) . Upon multiplication on CS phase with A = 3 
it leads to the following expression for the minimal spin 
component 



N 
i<j 



(47) 



notice the factor of Sim in the exponent which ensures 
the correct total density n. The kinetic energy of such 
state is factor of three higher than that ascending from 
Eq. (|4"Tj) . On the other hand, the correlation holes are 
wider, which may lead to a gain in the interaction energy 
(especially since the interactions due to higher spin com- 
ponents are essentially non-local). While at the moment 
we are not aware of a model where the fractional state is 
advantageous, it is worth pointing out that such a model 
is in principle possible. 

Notice that symmetric minimal spin component may 
be written with an arbitrary integer exponent p as a 

JV N 

-pirn J2 \ij\ 2 /2 ik J2 x j 



n 



JV 
i<j 



However only odd 



p's may be traced to CF construction. This seems to 
indicate that states with even p's yield larger average 
energy. This is indeed the case for p = 0, which is nothing 
but TRSB condensate with /ig a n. It is not clear at 
the moment how to demonstrate this statement for p = 
2,4 



E. Conclusions 

We have shown that the low-density phase of Rashba 
SO bosons may be described as the composite fermion 
state in the quantizing magnetic filed. Such a state is very 
different from both TRSB Bose condensate and SDW 
state discussed before. In particular its equation of state 
fi(n) oc n 3 / 2 leads to a different scaling of the kinetic 
energy and atomic cloud size with the number of parti- 
cles in shallow traps. It also implies different profile of 
the cloud density. The excitation spectrum is predicted 
to be gaped in the bulk of trap with the gapless chiral 
surface mode at its edge. For deeper traps we predict 
the phase separation between denser condensate phase 
in the middle of the trap and dilute CF phase at its edge 
with the first order density jump at the interface between 
them. 
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Appendix A 

Here we give details of calculations for the compo- 
nents of the two-particle wavefunction, ^-|'-|-(r, r) and 
*t4,(r,r) = * it (r,r) in Eq. ([26]) and derive Eqs. ([2?]> . 
(12811 of the main text. 



1. Calculation of ^^^(r, r) and the interaction 
form- factor fki,k 2 

From Eq. ([23]) we have 

$r tt (r,r)= (Al) 

where the fcrmionic part of the wave function is given by 
1 



*f (ri,r 2 ) 



V2V 



e ikiri+ik 2 r 2 _ gik 2 ri-Hk 2 r 2 



(A2) 



and anisotropic vectors f ^ and f' 2 in the Chcrn-Simons 
factor e lAarg ( ri_r2 ' arc defined as r'j = (aXj,y'j/a), j = 
1, 2, where a is the anisotropy parameter. 

To represent the integrand in Eq. (|A1[) in a convenient 
form, we convert it to the momentum space. Fourier 
transformation yields 



e iAargf = _ 2ni 



dk e aargk 



ikr 



(2tt) 2 k 2 



(A3) 



with k = (k x /a, ak y ). Fourier image of 7Z(r) is 7Z(k) 



/o -iar g k_ 
we obtain 



Substituting now Eqs. ([Mj) and (fM)) into ([ATI) . 



jgi(ki+k2)r r g iAargk 

*tt(r, r) = t= — / dk ■ 



k 2 



(A4) 



X I e _i ar s( k i+ k )- iar s( k 2-k) _ e -iarg(k 2 +k)-iarg(k!-k) 



In order to evaluate -Fkik 2 from (jA4|) . it is convenient to 
introduce the following notations 

e iAargk = k x /a + i\ak y = / z_ x 
|fc| \\2 
e - iarg(kl+ k) = ™i±* ^ = k + lk 
\wi + z\ 

p -i« g (k,-k) = _©2^* w 2 = k 2 , x + ik 2 , 



\w 2 - z\ ' 

Z — &x ~1~ I'ky • 



(A5) 



The complex variable z = k x /a + iak y can be conve- 
niently recast as z = |(a + 1/a) + | (1/a — a). Em- 
ploying integration variables z and z and using complex 
representation of vectors ki and k 2 , one obtains 



Ft 



kik 2 



% I dzdz 1 
V2. 



4ir \z\ 



- 1 



WIW2 



l\ z' 2 1 I z{wi-w 2 ) i 2 

y I wiiU2 W1W2 ' 



(u>l 



W1W2 



w 2 ) 



.(A6) 



In Eq. (|A6[) the short hand notation (wi « — >■ W2) stands 
for the same first expression in square brackets but with 
interchanged variables w± and w 2 . 

Upon introducing new dimcnsionless variables 



Wi - w 2 



y/WiW 2 



Eq. ([A6J yields 



kik 2 



WiW 2 



y/WiW 2 



dvdv 



(A7) 



\/2 V |u>lU>2| 7 4-7T 



^|P 2 -l + Z?f| 2 



(A8) 



(r -r) , 



where P = ^(a + 1/a) + ^(1/a — a) and parameter t can 
be rewritten in terms of original momenta ki and k 2 as 
follows 



(arg(ki)-arg(k 2 ))/2 



2 -i(arg(ki)-arg(k 2 ))/2 
fcl 



(A9) 



The integral on the right hand side of (|A8[) is convergent, 
therefore one can expand the integrand in f and then 
perform the integration. Keeping only linear in f terms 
in this expansion one finds 



^ki,k 2 = ic 
where 

c\(a) -- 



WiW 2 \ 



,-iarg(k 2 ) 



(A10) 




1 (|H 2 -i)(^ 2 -i)(^ + ^) 



The functions |c±(a)| are plotted in Fig. [7] 



(All) 



2. Calculation of ^-^(r, r) and the interaction 
form- factor G^ x ^ 2 

From Eq. ([23]) one finds 

tt n (r,r) = ± J dri^(r-ri)e iAaJg ^- ? )^(ri,r), 



(A12) 
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FIG. 7. (Color online) Interaction coefficients |ca(q)| and 
|dx(a)|. Full lines correspond to A = — 1 and dotted lines - 
to A = +1 chirality. 



where \I r F(r' 1 ,r) is defined in Eq. (|A2[) . while spin-rising 
operator 1Z is given by Eq. ([Tj) . After substituting these 
expressions into (|A12[) we obtain 



*t4.(r>r) 



i(ki+k 2 )r 



2V 



G 



kik 2 - 



v/2tt 



dr>- 



D iAarg(r-r')-iarg(r-r') 



(A13) 



r — r' 



l\2 



x |e ' ?;k 2( r - r ')_ e -' ik i( r - r ' 



which shows that Gkik 2 is a difference of two functions 
G kl k 2 = Iki - 4 2 - Introducing (r - r') x + i(r - r') y = 



re^, and j3 = log(a), we will have for A = +1 chirality 



D iarg(r-f 



cosh[/3]e^+sinh[/3]e-^ 
| cosh[/3]e^ + sinh[/3]e"^| 

cosh[/3]e^ +sinh[/3]e-^ 
^cosh[2/3](l + tanh[2 / 8] cos[20]) ' 



(A14) 



Importantly, the form-factor corresponding to the A = 
— 1 chirality can be obtained by interchanging cosh[/3] 
with sinh[/3] in the numerator of Eq. (| A14|) . Then we 
obtain 



1 



7r^2cosh[2/3] 

dr , (cosh[/31 +sinh[/3]e- 2 ^)e- 8fcrcos ^- arg(k )l 

— dq> 

r 



(A15) 



y/i + tanh[2/3] cos[20] 



We note that the integral Ik has a logarithmical divergent 
contribution coming from r — > 0, but it cancels out in the 
difference 1^ — /k 2 and therefore in Gk!k 2 - In order to 
find the coefficient d\ (a) in Eq. (|2"5)) we need to evaluate 
integral (|A15[) in linear over arg(k) <C 27r approximation. 
Upon expanding integrand over arg(k) and integrating 
over r we arrive to the following expression 



d\(a) 



d<t>- 



1 - cos [20] 



7r x /2cosh[2/3] Jo x/l + tanh[2/?] cos[2</>] 

(A16) 



where s + = sinh[/3] and s_ = cosh[/3]. The function 
|d±(a)| is plotted vs In a for both chiralities A = ±1 in 
Fig.El 
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